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SUMMARY 


Solutions are presented for the lift and pitching moment on a two- 
dimensional flat plate undergoing a step variation in pitch or vertical 
displacement. Resulting indicial lift and pitching-moment curves are 
given for free— stream Mach numbers equal to 0, 0.8, 1.0, 1.2, and 2.0. 
Considerable use is made of the analogy between the boundary values for 
a two-dimensional wing in unsteady motion and those for a three- 
dimensional, lifting s\irface in steady motion. The incompressible, 
unsteady case, for which Wagner’s classical treatment already exists, is 
shown to be analogous to a problem in slender wing theory. 


INTRODUCTION 


The response in lift for a two-dimensional airfoil that starts sud- 
denly from rest and moves forward at a constant velocity in an income 
pressible fluid was studied originally by Wagner (reference 1). It is 
well known that such a response, referred to as the indicial lift due to 
angle of attack, can, by suitable superposition, be used to find the lift 
on a wing undergoing an arbitrary variation of angle of attack with time. 
In particular, Garrick (reference 2) showed that Wagner’s results were 
consistent with those obtained by Theodorsen (reference 3) for a harmoni- 
cally oscillating wing. 

Subsequently, aerodynamicists have become interested in the effect 
of compressibility on the unsteady motion of a wing. Most, of this inter- 
est, however, has been focused on the problem of wing flutter which 
involves harmonic motion. Hence, the analysis of the compressibility 
effects has usually commenced with the immediate assumption that the 
motion is harmonic. Such studies which-, /iii general, can be considered as 
an extension of the work of Theodorsen, have been particularly successful 
in the case of supersonic flow (e.g. , reference 4). In subsonic studies, 
however, these methods have proven to be . s omfewhat cumbersome. 
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The purpose of this report is to calculate the effects of compress- 
ibility on the indicial responses. At supersonic speeds (e.g., refer- 
ence 5) some of these indicial curves have been obtained. The material 
presented here, however, will provide a complete set of indicial 
responses supersonic speeds and, what is 

more important, will also provide these same responses for subsonic 
speeds; the particular case of a free— stream Mach number equal to 0.8 
having been chosen for detailed consideration. Results are aJLso pre- 
sented for a wing traveling at the speed of sound. 

The analysis required to formulate these indicial curves is based 
on the two-dimensional wave equation having boimdary conditions speci- 
fied over the plane in which the wing moves. The particular boundary- 
value problems involved are, in fact, completely analogous to certain 
three-dimensional, steady— state, supersonic, lifting-surface problems. 

The latter field has received much attention diiring the last few years 
and analytical techniques have been developed by means of which many 
supersonic lifting-surface problems have been solved. Because of the 
analogy, these techniques can be applied to the field of two-dimensional, 
unsteady flow. The analogy even extends to the point where a two- 
dimensional wing moving unsteadily in an incompressible, field has for 
its analog a three-dimensional, lifting surface moving steadily in a 
compressible medium at, the speed of sound. By means of this correspond- 
ence, Wagner*s classical result can be rederived as a problem in slender- 
wing theory. 

The part of the analysis pertaining to the response of an unsteady 
wing traveling at a subsonic speed is lengthy and somewhat tedious 
re^rdless of the method of approach. With the use of indicial functions, 
however, the calculations are reasonably straightforward, especially for 
Mach numbers around 0.8 to 1.0. Further, the use of indicial functions 
sheds considerable light on the manner in which ^fcich number variations 
affect the section aerodynamic characteristics. 

LIST OF IMPORTANT SYMBOIS 


a. 


'm 


a. 


velocity of sound in the free stream 
chord of wing 
section lift coefficient 




section moment coefficient 


m 




indicial section lift coefficient due to angle-of-attack 
change (without pitching) 
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Cl ’ indicial section lift coefficient due to pitching on a wing 

'I rotating about its leading edge 

Cjj, • indicial section pitchingnmoiiient coefficient due to angle-of— 

^ attack change (without pitching) measured about the leading 

edge and considered positive when the trailing edge is 
forced down 


Cjjj « indicial section pitching-moment coefficient due to pitching 

^ measured about the leading edge for a wing rotating about 

its leading edge; considered positive when the trailing 
edge is forced down 



indicial section pitching-moment coefficient due to pitching 
measured about the leading edge for a wing rotating about 
its 3A-chord point 


I 

m 

Mo 

P 

Ap 

q 

<10 

r,s 

t» 


section lift 
section moment 
fli^t Mach number 
static pressure 


PZ - Pu 

v_ 


dsmamic pressure f ^o^o^^ 
oblique coordinates f r = s 


s(l+Mo) - cj2 


/t 


t+x \ 


1-4!, 


time 


t 

t 


apt* ^ t 
c c 


u,v,w 


perturbation velocity components in the x,y,z directions, 
respectively 


Au 


"u - ^Z 
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Vq flight velocity 

Xfj,z Cartesian coordinates 

o c 

a angle of attack in radians 

s J\ I 

r circulation strength 

5 function, normalized with respect to t’, thus 

&(t*)dt» = 1 


6(f) Dirac 



.(0) ^(f) 

9 angle of pitch in radians 

Pq density in undisturbed air 

cp perturbation velocity potential 


Subscripts 


u 

upper 

side 

of 

z=0 

plane 

2 

lover 

side 

of 

z=0 

plane 


1 variable of integration 

6 part of a response multiplied by a 6 function 


METHOD 


A basic linearized form of the partial differential equation which 
governs the, flow field surrounding a thin wing moving through the air 
can be written in terms of the pertvirbati on velocity potential, cp , as 


^xx ^yy '^zz 


_ 1 


^f f 


( 1 ) 
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where 


x,y,z 

spatial coordinates 

t» ' 

true time 

^0 

speed of sound 


Equation (l) is applicable wben the fluid at infinity is at rest with 
respect to the x,y,z coordinate system and the wing or body traces 
certain space curves for which time is the parameter; moreover the per- 
turbation velocity components u=<p^, ■v=(fy, w=(p2 must be small relative 
to the fli^t velocity Vq of the wing. 

In the following analysis a flat— plate wing moving in the z=0 
plane away from the origin along the negative x axis will be studied. 
Since the variations in the flow along the span will be neglected and 
the slope of the stream lines at the wing surface will be taken as the 
ratio of w^, the vertical induced velocity in the z=0 plane, ^ to Vq, 
the constant forward velocity of the wing, the analysis is consistent 
with the assvnnptions of two-dimensional thin-airfoil theory. Hence, for 
the incompressible case, the partial differential equation reduces to 

*Pxx ^zz ~ ® (2) 

since the speed of sound is very large with respect to the velocity of 
the wing, and for the con5>ressible case it reduces to 

cp +9 =9.. (^) 

XX zz ^tt 'J-' 

where t equals aot ’ . 

The boundary conditions reduce to the specification of w^ over 
that portion of the xy plane occupied by the wing as time progresses. 

The solution in terms of 9 can be converted into the loading 
coefficient by means of the equation 

^ = JL iSL = JL_ (4) 

qo Vq 2 at* VqMo at ^ ^ 

The results for section lift and pitching moment can then be evaluated 
by appropriate integrations. 

^he subscript u indicates the value of w as z approaches zero 
from above. Since the discussion is limited to flat— plate surfaces 
without thickness, this is the same as the value of w obtained 
as z approaches zero from below, and thus w^ simply refers to the 
value of w in the z=0 plane. 
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Boundary Conditions for 




(o) 

namely: 


the Indie ial Functions 

The . fundamental boundary— 
value problem to be considered in 
this report is the one generating 
the so-called indicial curves for 
loading, lift, and pitching moment. 
By definition, an indicial function 
is the response to a disturbance 
which is applied abruptly at time 
eqmls zero and is held constant 
thereafter; that is, a disturbance 
given by a step function. For 
example, if the angle of attack 
of a wing varies with time as 
shown in sketch (a), the resulting 
lift coefficient, also shown in 
sketch (a), is designated as the 
indicial lift coefficient due to 
angle of attack. Four such indi— 
cial functions will be evaluated. 


c^ the indicial section lift coefficient due to angle— of— attack 
^ change (without pitching) 



the indicial section pitching-moment coefficient due to angle— of— 
attack change (without pitching) measured about the leading 
edge and considered positive when the trailing edge is forced 
down 



the indicial section lift coefficient due to pitching on a wing 
rotating about its leading edge 


c • the indicial section pitching-moment coefficient due to pitching 
^ measured about the leading edge for a wing rotating about its 

leading edge and considered positive when the trailing edge is 
forced down 


The equations which transform these functions to those for a wing 
pitching about a point a distance ac back from the leading edge and 
having its moment center a distance be back from the leading edge are 
simply 
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‘'a ‘'a 


c„ = c„ • + t) c- 


■'w, 


h 


''a 


Cj =c, ac. 


= c„ ’ + b — a Cjj^* — ab cj 




a 


The boimdary condition which applies to the indicial functions due 
to angle of attack is simply that 


w^ Vq(x 


(5) 


over a certain planar area in xzt 
space representing the area occupied 
by the wing, and elsewhere pressure 
is continuous. For a wing of chord 
length c, this area, shown shaded 
in sketch (b), is boimded by the 
traces of the leading and trailing 
edges, X = — Vgt* and x = c— Ypt’, 
and the line along which the motion 
started (i. e., the x axis). 

In. the case of the indicial 
function due to pitch for a wing 
rotating about its leading edge, the 
boundary condition requires that the 
upwash be given by the expression 

Wu = -(x + Vot') 9 (6) 


•C— H 





over the same region in the xt plane as for the angle of attack case 
and, a^in, that pressure be continuous elsewhere. The angle of pitch, 0, 
is taken as positive when the trailing edge is lower than the leading 
edge, and 0 is the time derivative, d0/dt*, positive, when the trailing 
edge is falling with reference to the leading edge. 
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Direction of wing motion 


Flight path of leading 



(cj 


0 



The difference between B and a 
is illustrated in sketch (c). The angle 
of attack a is the angle between the 
flat, wing s-urface and the tangent to the 
flight path of the leading edge. The 
angle 9 is the angle between the flat 
wing surface and the horizontal. 

Sketch (d) shows a wing undergoing a 
sinusoidal angle-of-attack variation 
with a zero angle of pitch throughout. 

For convenience and in order to distin- 
guish from the pitching wing, such a 
wing will be referred to as a sinking 
wing. Sketch (e) shows a wing undergo- 
ing a sinusoidal angle— of— pitch varia- 
tion taken about the leading edge at a 
constant (zero) angle of attack. 

The variable q to which the lift 
and pitching-^aoment coefficients of the 
pitching wing are referred is equal to 

a dimensionless form of expression 
for rate of pitch. 


Adaptation of the Indicial 
Functions to Maneuvers 


fdj 



When the wing is undergoing a maneu- 
ver in which a and 0 make small but 
arbitrary deviations from zero, the bound- 
ary condition in the region of the xt 
plane traversed by the wing becomes 


= -Vq a - (x +,7Qt*) 0 (7) 


Since the theory is linear, the lift and moment on the wing are calcu- 
lated from the indicial functions by the principle of superposition. 
For example, the lift and pitching-moment coefficients developed by an 
arbitrary variation of a and 0 are 


'm 


o ^ 

t» 

[°ma’ (t'-ti’) a(ti')dtx* + Cn,^* 
o 


(t‘-tl‘) 


q(ti‘)<iti 


! 

t 

;(8) 


(t»-ti») q(ti')dti» 


J 
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Eqmtion (8) is derived for an angle-of-attack variation measTired 
with respect to the flight path of the leading edge and the primes on the 
derivatives indicate that the wing is pitching about and the moments are 
measTored with respect to the wing leading edge. If the angle of attack 
is measTired with respect to the fli^t path of some other point on the 
wing (as, for example, the center of gravity), then it can he shown for 
small deflections that 


a = ai — 


ci9: 


0 = 01 


where ai and 0i are the angles measured with respect to the new 
flight path and c5 is the distance from the leading edge to the new 
reference point. (See sketch (c).) The total value of the section lift 
and pitching-moment coefficients for an arbitrary variation of ai 
and 01 would be 


cz = 





dt» J 


o 



0'l(tl’ ) 


0 'l(ti*) 



- tqi(ti’) 


+ C2g*(t‘-ti»)qi(ti‘)J> dti« 

i+ •(t’-ti»)qi(ti»)l iti* 
J ^ J 


Adaptation of the -Indicial Functions to 
Flutter Derivatives 


The notation adopted in this report does not coincide with the 
usual notation used in the study of fluttering wings; however, it should 
be sufficient to relate the present notation to the one used in refer- 
ence 4. Thus, the values of the terms Li, L 2 , Lg', L^* and Mi’, M 2 *, 
Mg*, M 4 * used in reference 4 can be evalmted as follows: 


ifa)t* H r / ^ iuti* 

— 'icjt * t * 

Mi*+iM2* = F c *(t*— ti’)e^ dti* 

2k t*->» dt*J 


) ( 9 ) 


o 
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— iwt’ 


^3 ^^'^4 


t * 

[c!a('^'-‘i')H<- c,^'(t'-t,*) 


e dti' 


>(9) 


w lim d f 


J^7J ’i 




THE SOIOTION FOR THE IWDICIAL FUNCTIONS 


The boundary conditions which represent a wing undergoing a sudden 
jump in pitching velocity or angle of attack have been discussed and 
presented in equations (5) and (6). Their application to wings flying 
at several Mach numbers will now be presented. The incompressible case 
will first be discussed, and then the cases in which M© equals 0.8, 
1.0, 1.2, and 2.0 will be considered. 


Incompressible Case, Mq = 0 


The solution for the incompressible case is applicable when the 
forward speed of the wing is small compared with the speed of sound so 
that the ratio can be neglected in comparison to tinity. The 

basic partial differential equation governing the flow field was presented 
in equation (2) as 


*^zz ® 

subject to the boundary conditions discussed, and the equation for the 
loading coefficient was given in eqiiation (4) as 

_ ^ ^ 

° Vo^ 8t« 

This boundary— value problem corresponds precisely to that which is 
studied in three-dimensional, steady— state, lifting— surface theory under 
the classification "slender-^ing theory." (See references 6 and 7.) 
Such an analogy is useful since well-established concepts in one theory 
can be immediately carried over into the other, (it should be re- 
emphasized, perhaps, that the subsequent treatment of the incompressible 
case is not intended to be an improvement on Wagner’s original deriva- 
tion but rather it is a rederivation along lines that will be used later 
in the analysis of the compressible cases. ) 
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The initial piilse .— The first analog 
vith s lender-wing theory which will he used 
concerns the initial pulses that occur in the 
values of lift and pitching moment. It is a 
well-known result (reference 6 ) that the 
total lift, as given by slender-wing theory, 
on the wing shown in sketch (f) is a func- 
tion only of the maximum span and the value 
of w^j along the section of maximum span 
(section AA). It is, therefore, independent 
of the wing twist and leading-edge shape 
ahead of section AA. This concept has been 
extended in slendej>^ing theory to the extreme 
case shown in sketch (g) of a rectangiilar 
wing. The lift on such a wing is concen- 
trated entirely along the leading edge and is 
a function only of the span of that edge and 
the value of w^ there. By the analogy 
existing between the two theories, therefore, 
it is evident that the solution to the indi— 
cial problems in two-dimensional, incompress- 
ible, unsteady flow (sketch (h)) will contain 
a pulse at t’= 0 . 

The evaluation of this pulse will be 
treated briefly. A solution to equa- 
tion (2) for the vertical induced velocity 
in the z =0 plane can be written in terms 
of the jump in u across the z =0 plane 
(see reference j), thus, for the shaded 
area in sketch (i) this is 

»u<>=) ■ - S ^ (10) 

The general inversion of eqmtion (lO) 
can be written 


y 



Pulse in lift at O. 




Au(x) 


?r V (x+a ) (b-x ) (x+a ) (b-x ) 



w,i(xi) rz : 

hj (xi+a)(b— xi) dxi (11) 

X— Xi 


where 



-a 


Au dx 
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In the present case A is zero, since A<p is zero at x=— a and x=h, 
and an integration of both sides of eqmtion (ll) with respect to x 
between the limits b and x gives 


Acp(tSx) 





J (1)-xi)(a«) *■ 

Ji x-xi)(a+b) 


Adoption of the notation 

v(t*,xi) + constant = w(t*,xi) dxi 
and integration by parts leads to the equation^ 


(T>-x)(a+x) 

Jt 



'^(tSxi) dxi 

(x— Xi )iy(b-xx)(a+xi) 


( 12 ) 


The loading can now be determined by using equation (4). If the 
shaded area in sketch (i) is allowed to vanish, all the loading accimiu— 
lates along the x axis in the region 0 £ x < c. Therefore, the inte- 
gral of the loading with respect to- t* over the shaded portion must be 
considered. The final result for the pulse loading (^/qo) t'=0 

can be expressed in terms of the 6 function (see list of symbols) as 


V^o/ J (x-: 


v(0,xi) dxi 


•xi C-Xi) Xi 


(13) 


The boundary conditions for the sinking and pitching wing given by 
equations (5) and (6), when inserted into equation (l3), yield 




( 

' q 6 / 


( 14 ) 


o 5 


The constant gives zero when placed in equation (12 ), provided 
b > X > — a. 
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After integration, the pulse values for lift and pitching moment 
may also he obtained. Hence, 



where the primes indicate that the wings are pitching about and the 
moments are measured about the leading edge. These expressions may be 
inserted in equations (8) and (9) and, since the integrand becomes zero 
everywhere except at the point ti*=t’, the expressions for the lift and 
pitching^noment coefficient developed by an arbitrary variation of a 
and G with time are 



c 


I 

m 


itc ^ _ 9_^ 
Wq 6h V 


The variation for t* > 0.— The integral equation (lO) is still per- 
fectly valid when applied to the flow field for t’ > 0. It is convenient 
to rewrite the equation in this case, however, so that the effects of the 
vorticity on the wing and in the wake are separated. Thus, 


'^u = 


r°“^o‘t' Au(tSxi) 
2it J X— Xi 

-V t» 
o 


2k 


Au*(xi) 


X— Xi 


dxi 


(16) 


c-V^f 
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(j) 


where Au*(xi) is the value of Au in the 
starting vortex wake. It is independent 
of t* since its value at all points along 
the line ab in sketch (j) is the same as 
at the point a. 

A reduction of equation (l6) can be 
obtained for the case of the sinking wing, 
where w^= by using the inversion 

given by equation (ll). Thus, 


JtAu(t’,x) J (x+V-t * ) (c— V ± X ) = A-Jt(2x+2V t’-c) V a + 


J Au*(xi) - 


c-Vot* 


r y(xi+Vot»)(x-^+Vot«) 

(x-xi) 


dxi 


Since A is given by the relation 


c-V^t‘ 

A = / Z!ai(t’,xx) dxi = — 


Au*(xi) dxi 


-V 


c-Vot* 


it follows that 


Voa(2x+2Vot*-c) 

Au(t«,x) = 2 L_ + 


J (x+T^t * ) (c-V^gt »-x ) Jt y (x+V^t • ) ( C-Y^t *-x ) 


pC y(xi+Vot»)(xi^+Vot«) dxi 

J 

c-V t» 
o 


(17) 


According to the Kutta condition Au(t’,x) vanishes as x approaches 
c— V^t*; the integral eqmtion for Au*(xi) thus becomes 


Vrft = 


-i/ 


Au^(xx) /jSilZot: 


c-Vot» 


Xi^+Vgt* 


dxi 


(18) 
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The section lift and pitching moment can be derived in terms of 
Au*(xi) in the following manner. By definition, the section lift I 
is 


I 



c-V«t* 


Po^o 


-V* 


^dx = 
'lo 


Po 


pc— V 

-v.t» 


c-Vof 


^(p 

at‘ 


dx 


(19) 


Since the value of AT is zero at the leading edge and at the trailing 
edge is equal to the total circulation F, two alternative forms for the 
lift can be written 


^ = Po^o^ + ^o~^J Acp(t‘,x) dx (20) 

-V’ 

and 

c-Vot» 

l = xMf,x)dx (21) 


By substituting equation (17) into (21) and integrating, it can be shown 
that 


I = J' 


c 


Au*(xi) 


c-Vot* 


y (xi+Vot ' ) (xi-c+Vot ‘ ) 


dxi 


( 22 ) 


Since the value of Au*(xi) has 
been determined by Wagner, equa- 
tion (22) represents a solution 
for the section lift. A plot 
of its variation in coefficient 
form is shown in sketch (k). 
Initially there is the pulse 
having an intensity defined by 
equation (15). After the pulse 
at t’=0, the value of the sec- 
tion lift coefficient starts at 
one-half its asymptotic value. 

It then increases, slowly 
approaching its asjTnptote of 2n. 



fkj 
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By definition, the pitching moment can he written 


,c-V_t» 


m = 




-V t« 
o 


(23) 


,c-V _t» 


-Pof ° (Vot’+x).^dx 


where the moment is taken about the leading edge and is positive when 
the trailing edge is forced down. A development, similar to the one 
given for the lift, gives 



(24) 


This result, that for t >0 the indicial center of pressure remains con- 
stant at the quarter chord throu^out the motion, is classical. 


If the boundary condition for a pitching wing, w^ = — (x+Vot*) 9, 
is substituted into equation (l6) and the inversion given by equation (ll) 
is again used, it can be shown in the same manner used in the' derivation 
of equation (l8) that for x = c— ?ot’ the relation 




_i r 


jt J 


Au* (xi ) 


c-V t * 
o 


Xi +Vpt* 

Xi— c+V Qt * 


dX] 


applies. This integral equation applies to a wing pitching about its 
leading edge. If, instead, the wing is pitching about the three- 
quarter— chord position, an essential simplification is achieved. In this 
latter case, downwash is given by the expression 


w^ = -0(x+Vot»-| c) (25) 


and the resulting integral equation becomes 


0 = 


AUg*(xi) 


c-Vpt* 


Xi-c+V_t» 


: dxi 


(26) 
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where Au 3 *(xi) represents the vorticity in the wake following such a 
motion. The solution to equation (26) is simply 


Aua* (xi) = 0 


( 2 ?) 


From equation ( 27 ) it follows that the total indicial lift for t ’ > 0 
on a wing pitching about the three-quarter chord point is zero, and that 
the wing wake is free of vorticity. Further, it can be shown that the 
total' indicial pitching moment (still measured about the leading edge) 
is 


m 


pJtc3 

■ 15 " 


Vo 9 


( 28 ) 


The transfer of equations (27) and (28) back to the case in which 
the wing is pitching about its leading edge can be readily accomplished 
by means of the boundary conditions shown in sketch (z). Hence, if 
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By means of equations (24) and (29) the expressions for the three 
indicial functions, , 02^', and Cn^’ can all he written in terms of 
the indicial lift function for t>0. Hence, 


c ’ = 2- c 

^q F 

‘"V " “ ll ■ I 


> (30) 


The variations of the four indicial functions are shown in figure 1. 
For values of Vot*/c larger than those shown in the figure the approxi- 
mate equation suggested in reference 2 can be used, namely. 


^ 2+(VotVcT] 

This alternative result has, according to reference 2, an error of 2 per- 
cent or less for the entire range of time from 0+ to infinity. 


C7 = 2n 
‘a 


Subsonic Case, Mq = 0.8 


When the Mach number is no longer small, the analysis in the pre- 
ceding section must be modified. As an example of this modification, the 
calculations for a wing traveling at a Mach number equal to 0.8 will be 
carried out in detail. 

Equation (3) presents the basic partial differential equation of 
the flow field, and equation (4) gives the expression for the loading 
coefficient. The analogy which existed in the incon5)ressible case 
between the theory for the unsteady, two-dimensional wing and slender— 
wing theory exists in this case between the theory for the unsteady, 
two-dimensional wing and the theory for a steady— state , three-dimensional 
wing traveling at a supersonic speed. Thus, in the three-dimensional, 
steady— state case the partial differential equation governing the flow is 

9yy + 4*zz ~ 3 ^xx ( 31 ) 

and the equation for the loading coefficient is 
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Ap _ 4 

q, Vo ^ 


(32) 


The boxmdary conditions are in both cases that cpz is given over a 
portion of the z=0 plane. It is evident by a comparison of equa- 
tions (3) and (31) and equations (4) and (32) that results from the 
three-dimensional, supersonic, steady-state study (hereinafter referred 
to as the steady— state case) can be transferred to the two-dimensional, 
unsteady study (hereinafter referred to as the unsteady case) simply by 
replacing x, y, and 3 in the former case by t, x, and 1, respectively, 
and by dividing the result for the 
loading coefficient by Mq. 


The analog to the boundary 
condition for the problem of 
finding the indicial loading on a 
tvo dimensional wing flying at a 
tfech number equal to 0.8 
(sketch (m)) is the boundary con- 
dition for the problem of finding 
the loading on a constant— chord, 
swept— forward wing tip with a 
subsonic trailing edge such as 
that shown in sketch (n). The 
Mach cones in the steady-state 
case, traces of which are shown 
as dotted lines in sketch (n), 
become, in the unsteady— state 
analog, the locus of the sound 
waves which started at t=0 from 
the leading and trailing edges of 
the two-dimensional wing 
(sketch (m)). Finally, the 
analog in the steady-state field 
of the unsteady wing would be a 
flat plate for the unsteady 
sinking wing and a plate with a 
linear variation of twist for 
the unsteady pitching wing. 

A detailed analysis of an 
unsteady, two-dimensional wing 
flying at a Mach number equal to 
0.8 will be presented. Just as 
in the section on inconpressible 
flow, the study will be divided 
into two parts. In cases for 
which Mq^O, however, the indi— 
cial functions contain no pulse 
at t=0. Hence, the first part 



(m) 



(n) 
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of the Mq= 0,8 study will he concerned with the behavior of the indie ial 
functions in an interval for which t is small but finite (actually 
0 ^ t ^ 5c) and the second part, with their asyngjtotic behavior. 


The early stage.— The analog vhich exists hetveen the steady-^tate 
and unsteady cases may be utilized to great advantage since the large 
number of special methods and techniques developed for the solution of 
problems in the former case may be applied to the solution of the 
analogous problems in the latter field. In this manner an exact solu- 
tion for the loading over the five regions shown in sketch (o) may be 
obtained for both the sinking etnd the pitching wing in an unsteady flow 

field by the use of methods such 
as those presented in reference 8. 
Solutions for larger values of 
time could also be obtained^ but 
the labor involved in calculating 
such cases becomes prohibitive and, 
as will be shown later, approx- 
imate methods can be developed 
which extend the solutions for the 
indie ial lift and pitching-moment 
curves to their asymptotic values. 

The analysis used to calculate 
(o) the loading in terms of Xq = x/c 

and to = t/c over the five 

regions shown in sketch (o) on both the sinking and pitching wings is 
outlined in appendix A. Plots of the loading on a sinking wing are 
shown in sketch (p) and, in more detail, together with the loading on a 




pitching wing, in figure 2. At to=0 the 
loading is constant,® and as time increases 
the loading-coefficient ctirve approaches the 
familiar two-dimensional, steady— state shape 
given, for the sinking wing, by the equation 

^ = (33 

q© P V X 


(p) 


Q 

The result that the initial shape of the load distribution is the same 
as the shape of the given curve for w^ also applies to all three- 
dimensional wings of arbitrary plan forms traveling at subsonic or 
supersonic Mach numbers. 
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where x is the dlstejice from the lesding edge which is at the point 
x=0. Sketch (q) shows that the distribution at to = t/c = l/(l-Wo) = 5 
is already essentially the same distribution as that obtained at to=» 
(i.e., the agreement is good with the curve produced by multiplying the 
right side of eqviation (33) by a constant factor).^ The use of this fact 
simplifies subsequent analysis concerning the asymptotic behavior of the 
indie ial lift curve. 



The indicial lift and pitching-^noment functions were also calculated 
(see appendix A) in the range 0 ^ to ^ l/(l-^lo)* Their variation in 
this interval is shown in sketch (r) for Mo = 0.8. It is evident that 
the calculations must be extended beyond the point to=5 since the 
asymptotic values are not 3 ret even closely approached. 

Before studying the nature of these curves for large values of to, 
however, it is xiseful to examine them with reference to the disciission 
in the previous section on incompressible flow. For example, it was 
pointed out that the indicial center of pressure on the sinking wing 
remained at the quarter— chord point for t’>0. It is, therefore, pei>- 
tinent to consider the location of the center of pressiire on the sinking 
wing when Mo = 0.8. By means of the curves given for ci^ and Cn^* 
and by the relationship 


4 


A similar result was noted in the study of the load distribution on 
swept-back wings with subsonic leading edges (reference 9)* 
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the variation of (x/c)j,^p^ is easily evaluated ((x/c)g^p^ is the dis- 
tance between tiie leading edge and the 'center of pressure divided by the 
total wing chord). !This variation Is shown in sketch (s). It is 
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apparent that the center of pressvire is very close to the quarter— chord 
point for values of to>5- Hence, at M© = 0*8 the significant effect 
of compressibility on the location of the center of pressure is contained 
in the interval 0 < t© < 5- This leads, immediately to the further use- 
ful conclusion that for values of to greater than 5 the value of the 
pitching— moment coefficient is given essentially by the equation 

Likewise, it is apparent from the discussion of the incompressible 
case that the indie ial functions for the pitching wing can also be 
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expressed in a more convenient form by shifting the axis of rotation 
from the leading edge to the three-quarter— chord point. Using the values 
of Ctt^* . Cj *, and Cm^' presented in sketch (r), the variation of 

^®Zq^ 3/4 calculated from the definitions given in 

equation (29). These curves are presented in sketches (t) and (u) and 
again it is apparent that at a Mach number equal to 0.8 the compressi- 
bility effects axe contained in the interval 0 < < 5* 

The later stage .— It follows from the preceding discussion that when 
to is large, the values of the indicial functions 

for compressible flow can also be expressed adequately in terms of c^ 
by equations similar to equations (30) which were derived for incompress- 
ible flow. Thvis for to > 5> one can write 

'1 

Cma' = - 

'^^q* = 3 Hjh I (34) 

cn^’ = - (3c^^/l6) - («/8p) 

The four indicial responses have all been shown to depend only on 

the value of the function Cj if to is large. It remains, therefore, 

oc 

to determine the asymptotic behavior of This can be accomplished 

in the following way. Consider the steady— state solution for the lift on 
a two-dimensional, flat lifting surface traveling at a subsonic Mach 
number. As was pointed out by Wieghardt (reference lO), if the lift on 
such a surface is represented by placing at the quarter— chord point a 
vortex which has the same circulation as that developed by the wing, the 
angle of attack measured at the three-quarter-chord point will be the 
same as that of the simulated lifting surface. Extending this concept to 
include the unsteady effects, an investigation will be made of the varia- 
tion with time of the vortex strength which will maintain a constant 
angle of attack at the three-quarter— chord station following an impulsive 
start at to=0. 

The analogovis problem in steady— state theoiy becomes one of finding 
the strength of the vortex system, shown in sketch (v), which gives a 
constant value of w along the line CD.^ 

Each vortex composing this system lies along the line AB, extending 
from the minus infinity toward the origin, and trails back parallel to 

^In the vicinity of the origin, of co\irse, this representation gives a 
poor approximation to the original boundary— value problem; hence, use 
of the results must be limited to the regions in which to is large. 
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(v) steady ~ state analog Unsteady case 


the xi axis to form the trailing vortex sheet. Note that, for con- 
venience, the origin of the axis system has been located at the quarter- 
chord point. The solution to such a problem in steady— state, lifting- 
line theory would result from the solution of the integral equation® 

(for a development using the notation adopted here, see reference 11); 


w 



ACp(yi) dxi 

[(x-xx)2 - P2(y_yj2]3/a 


where A9 is not a function of xi since the strength of a trailing 
vortex is, of course, constant. The area of integration t is the 
region within the forecone springing from the point x,y. 

If the above equation is transformed by means of the analogy to 
represent the solution of the unsteady problem (see sketch (w)). 


'The symbols ^ and ^ are used to indica 
be taken. Thus (see reference 11 or 12), 


are used to indicate that the finite part is to 


(x-y)‘ 


dx Ja (x~y) 


where G (x,y) is the indefinite integral of f(y)/(x— y)^. Further, 


r = _2 

dt (x-y)'^/^ 


fill_ 


a V x-y 
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3, X, and y are replaced by 1, and x, respectively; and Acp, the 
total jun^ in potential at a given section, is replaced by the circula- 
tion r. Hence. 

where as indicated in sketch (w), is the area in the forecone from 
the point P which lies always along the linfe x = (l/2)c — Mot. 
Integration with respect to ti reduces the last equation to 


w 


2« J Mo(x-t) 

1 + Mo 



which, by means of the substitution Xi/c = X 2/2 becomes along the line 

X = I - Mot 


P^O (A-O+dO-Xs) r(x2) 

CK^ To (|l2-X,Q+X2)^ 0-X2 )(^-o'H^1~^2) 


dX2 


(35) 


where ko = 2 Moto - ^lo^ l^o = Mo/(l+Mo), = 2 Mo/(l-^>4o^) > = l/(l+Mo), 

to = t/c, and where, of course, w is a constant equal to — Vo<x. 

A solution for r(x 2 ) in the integral equation (35) inay be obtained 
by expanding P in a series of the form 


r(x2) 


= r / ^2 ^ ^ r ^2 ^ 

P Lv (ao+X2) ^ V (ao+X2)® 


(36) 


Place equation ( 36 ) into (35) and expand in powers of lAo- There 
results the expression 

w _ T . Cl . C 2 Zn(lAo) . C 3 

_ “■ ..L ' ^ ^ o ' _ Q > • • • 

VqCL Xq 


in which 
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Cl = bi - 


Sq 

2 


+ 


32 


Hence, if ag and bi are chosen so that Ci is eq.ua! to zero, an 
expression for T will be obtained which represents the solution to the 
integral equation (35) correct to the first order in. l/to (i.e., 1/Xq) 
for large values of to. Further, if equation (36) is expanded in 
powers of l/xg^ there results 


r(x2) 


rtcaVp 

3 


1 ^) + 



which becomes, using the condition for Ci and relating xg and tp 
by the equation of the leading edge. 


r(tp) 


jtcaVp 


1 - 


2Mpto32 



(37) 


The relation between circulation and lift in incompressible flow 
has been derived and presented as equation (20). For compressible flow, 
this expression becomes 


N pl-Mptp 

I = PpVpr + ppap / Acp(to,Xo)dxo 

0^0 ^Mptp 

where Xp = x/c.- In order to obtain a complete expression for the sec- 
tion lift, it is obviously necessary to know the chordwise variation of 
A<P. The as 3 miptotic behavior of I can be calculated, however, by 
applying the results presented in sketch' (q). This result suggests that 
for large values of tp the value of A9(tp,Xp) used in the equation 
for section lift can be expressed by the product of AP( op,Xp) and 
r(tp)/r(»). In other words, for large values of tp the shape of the 
chordwise distribution of vorticity is the same as the two-dimensional, 
steady-state value. An indication of the accuracy of such an approxima- 
tion is shown in sketch (x) where the precise value of A9(tp,XpJ is 
coDq)ared with the approximation at tp=5. 

Since r(») = «aVpc/3, the substitution of 

A'Ktp,Xp) = Ap(oo,Xp) r(tp)/r(oo) 


^ - Po^oT ^ PoS-0 ^ 


/ 


gives for the lift 
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(x) 


Place equation ( 37 ) in the above equation and it is evident that, to the 
first order in l/t©, the section lift is given by PoVoT. Hence, it 
follows that for large to 


c 





1 

2MotoP2 


+ • 



(38a) 


There remains the problem of joining the above result for with 

the one derived in the preceding section and valid for 0 <to^ 5 « To 
accomplish this end the equation 



1 hi 

ho + 2MotoP2 (ho + 2MotoP2)2 


(38b) 


was used to express ci for the range 5 ^to. Obviously the value of 
hy equation ^38b) has, to the first order in l/to, the same 

asymptotic variation as that given by equation (38a) regardless of the 
values of the constants ho and hi. These constants can be chosen, 
therefore, so that both the magnitude and slope of the indlcial cj^ 
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ciirve given by equation ( 38 b) are continuous at to=5 with the exact 
ct curve obtained for 0 < to ^ 5 in the previous section. 

OG 


The curves for f ‘^niq*^ previously discussed, can 

be calculated from cj by the use of equations (34). In each of these 
three cases, the magnitude of hi was modified slightly so the resulting 
curve would be continuous at tg=5 with the exact results presented in 
the preceding section. The final expressions are 


2n 


“■a 




131 -g " 






■ p L 

17. 06 + to 

( 17 *06 + to)^ _ 

“It 

\ 1.736 

121.8 

■ 23 

17.06 + to 

(17.06 + to)2 

_ 3 « 

[l 

134.3 

23 

17.06 + to 

(17.06 + to)^ 


1.302 

90.53 


17.06 + to 

(17.06 + to)^ 


> (39) 


J 


The final indicial section lift and pitching-moment ciorves are 
shown in figure 1 plotted against the .parameter Vot*/c, the number of 
chord lengths traveled (Vot’/c = Moto) . Tabular results of the indi— 
cial cmrves are also presented in table I. ' 


Sonic Case, Mq = 1.0 


The general results, obtained in the preceding section and pre- 
sented in appendix A, for the indicial loading over the sinking and 
pitching wing may be extended to the sonic case. Firrthermore , the two 
intervals for which analytic results in a closed form were presented in 
appendix A now cover the complete time range since 0 ^ to ^ l/(l+Mo) 
becomes 0 ^ to ^ O .5 and l/(l+Mo) ^ to — l/(l-Mo) becomes 
0.5 — to — <». Hence, by an appropriate limiting process equations (A8), 
(A9), (AIO), and (All) become for 0 ^ to ^ O .5 

Cmix* ~ ^ 4o^ 

C2q» = 2 + to^ 

cmq’ = -(Vs) - (2/3)to® 


(40a) 



30 


MCA TN 2403 


and for 0.5 ^ t© ^ “ 

= (4/«) (2 ^/2to-l 






+ arc. cos 


-) 

« = -(2/n) J2t^ +(l- 

.(1 + 


°Zq* = (2/it) 


Cn^' = (4/3rt) 


^ v^2to-l +(l + ^ 

l4-to-3to^ 


arc cos 


arc cos 


to-1 


to-1 

to J 






)(40b) 


arc cos 


to-ll 

to J 


J 


Since the magnitude of the functions in equation (40b) grows indef- 
initely with increasing time, the assunq)tions of linear theory are even— 
ttially violated. However, for moderate values of to, these f\mctions 
have the same order of magnitude as similar indie ial curves for Mch 
numbers other than 1. These effects are Illustrated in figure 1. 


Supersonic Case, Mq = 1.2 and 2 

The method of obtaining solutions for the indicial functions at 
supersonic Mach numbers parallels the development presented for Mq = 0.8. 
The steady— state analog to the supersonic unsteady wing problem is a 
constant-chord wing tip with a supersonic trailing edge. See sketches (y) 
and ( z ) . It is well known that the problem of finding the loading over 
wing plan forms with all supersonic edges is one of tfie simplest in three- 
dimensional, lifting-surface theory. In fact, since the upper and lower 
surfaces are noninteracting, the solution is determined by integrating 
sources within the Mach forecone. The analysis for c^^ has already 
been carried out in reference 5« 
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The analysis iised to calculate the loading in terms of Xq=x/c 
to=t/c over the three regions shovn in sketch (aa.) is outlined in 
appendix B. An exan^le of the manner in which the loading varies with 
time over a sinking wing traveling at a Mach number equal to 1.2 is given 
in sketch (ab) and, in more detail, together with the loading over a 
pitching wing, in figure 3. 




(abJ 


The expressions for the indicial lift eind pitching-moment coeffi- 
cients are given analytically in appendix B, and plotted in figure 1. 

It cem be shown that the results given in appendix B reduce to the 
expressions given by equations (40) when Mq is allowed to approach one, 
so that there is no discontinuity in the theory in passing through the 
sonic range . 


COHCLUDIBG BEMAEKB 


The use of the analogy between the basic flow equations in steady- 
state, lifting-surface theory and in unsteady-state, airfoil theory has 
resulted in a method of calculating two-dimensional indicial functions 
throughout the subsonic and supersonic flight range. The results are, 
of course, subject to the restrictions of linearized, compressible— flow 
theory and, for example, the calculated responses given in figure 1 for 
sonic speeds must be considered as being outside the realm of validity 
within a few chord lengths of travel. la application to high-frequency 
oscillations, however, the initial portions of the indicial curves 
dominate the response characteristics of the airfoil and calcxilations 
near M© equal to one need not be invalid. 

In the supersonic range, the expressions in appendix B are calcu- 
lated for arbitrary Mq. It is apparent that no fixed Mach number effect 
can be used in transient responses except at high values of flight Mach 
number. In the subsonic range, however, the expressions apply only for 
the period of time 0 < to ^ l/(l-Mo). For values of to ^ l/(l-^do)^ 



32 


MCA TW 2403 


the method outlined in the report appears to he satisfactory for Ifech 
nxunbers equal to or greater than 0 . 8 . Preliminary calculations made at 
Mach numbers less than 0.8 indicate that it is necessary to extend the 
exact solutions past t© = l/(l-Mo). This extension is feasible if the 
cancellation techniques outlined in reference 13 are employed and the 
more difficult integrals are expanded into series form. 

Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 

Moffett Field, Calif., April 12, 1951. 



NACA TN 21^03 


33 


APEEKDIX A 

DETEEMHIATION OF SUBSONIC, INDICIAL, SECTION LIFT AND 
PITCHING-MOMENT CURVES 
TEE LOAD DISTRIBUTION 


The following resiolts for the indicial load distribution on sinking 
or pitching wings can be obtained in two ways. One of these methods will 
be outlined in the subsequent paragraphs. The other is outlined in 
references I 3 and li<- and is referred to as the lift— cancellation techni- 
que. The latter method has been used to check the load distributions 
originally obtained by the former so that an independent check of these 


results has been carried out. 

It was shown that the lifting- 
surface analog to the solution for 
load distribution over an unsteady, 
two-dimensional wing traveling at a 
constant subsonic speed involved the 
calculation of load distribution 
over a swept— forward wing tip with 
subsonic edges. Sketch (ac) indi- 
cates the geometry associated with 
the boimdeiry conditions, and the 
solutions are calculated for the 
various regions shown. 



(ac) 


In the notation of the imsteady 

prohlem the expression for the velocity potential can be written 


1 P P (wu) dti dxi 


(Al) 


where t is the area on the wing plan form included in the Mach fore— 
cone from the point (t,x). Equation (Al) is applicable only for cases 
in which w^ is known at all points within the forecone, as is the case 
when the edge of the wing within the forecone is everywhere supersonic 
(i.e., region 1 in sketch (ac)). However, Evvard (reference 8) has 
extended the solution provided by equation (Al) to include cases such as 
shown in sketches (ad) and (ae) in which the forecone intersects a sub- 
sonic edge and includes a region of imknown upwash. As was pointed out 
in reference 8, equation (Al) applies in these instances if the area of 
integration t is limited to the shaded regions shown in the sketches. 
It is apparent, therefore, that the potential (and thus the loading) 
over a sinking or pitching' wing can readily be determined for regions 1, 
2, and 4 in sketch (ac). 
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Points in regions 3 and 5 have forecones which intersect two sub- 
sonic edges, and the method just discussed can no longer be directly- 
applied. In reference 8, however, a method was given of evaluating the 
upwash in the region between the Mach cone from the apex and the leading 
edge (region 6 in sketch (af)). Thus, the plan form has become, effec— 
ti-vely, one such as shown in sketch (ag) in’ which only one edge is sub- 
sonic. This reduces the problem of finding the potential in these 
regions to the same problem as was involved in region 4. The analysis 
used in finding the loading over the various regions will now be 
considered. 



(af) ( ag) 

First, introduce a new coordinate system in which the lines x = — t 
and X = t are taken as the r and s axes, respectively, (see 
sketch (ab)). (This amounts to a rotation of the original axial system 
through an angle of 45°.) The transformations relating the r,s to the 
x,y system are 
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1 / N 1 / N 

r = — (t-x) t = — (r+s) 

a/2 a/2 


s = (t+x) X = (s-r) 

V2 


In the new coordinate system equa- 
tion (Al) is written 



(wn)dridsi 

'(r-ri)(s-sx) 


The vertical induced velocity 
wu over the wing plan form is given 
in equations (5) and (6) for the sink- 
ing wing and pitching wing, respec- 
tively. The method developed in 
reference 8 was used to obtain the 



(Oh) 


value of Wu over the area between the lines s=0 and r=s(l+Mo)/(l— M q) 
(region 6 in sketch (af). The results for the sinking and pitching wing 


are, respectively. 



In terms of the r,s coordinate system, the expression for the 
loading can be written 


8-0 ^o^ \Sr Ss J 


(A4) 


If the subsonic trailing edge is not included in the forecone from the 
point at which the loading is being determined, the expression for cp 
may be substituted directly into equation (a 4). However, if the trailing 
edge is included in the forecone, and if the Kutta condition is to apply 
along such an edge, it can be shown that the equation for the loading 
coefficient assumes the form 


\ 



Hid Hid Hi d 
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^ _ 2 r dri f _i_ + ^ r (vu)<^si 

<lo itVoMo Jy J r-n N JJ s-si 


(A5) 


Sinking Wing 


The preceding method can be applied to the sinking wing to obtain 
the following integral .relationship for the loading over the various 


regions . The subscripts on 
particular equation applies. 


<lo / 


indicate the region for which the 


VWi 


2 ^ a ^ 




s J 




ds 


r-Ti 


■riy 


SHSi 


/ 

I 


^ol 


A. A 

itMo V Sr Ss 


r(l-Mo) a/ s— s 

(1+Mo) 


dsi r 

S—Sx J — 


dri 


-si V r— ri 


C^3 


2 J j + A 

a/ r-ri \ 


dS] 


ri(l-Mo) a/" S-^1 

(1+Mo) 


Jo J r-ri \^ri Ss/ Jo 


ri(l-^) 

(1+Mo) 


|dsi 




ri-si 


2si I 2s 1 

arc tan ^-r 


)^(ri+si) 


(ri“ Si)~Mo(ri+Si) 


a/ S-Sl 


° dri / a t ^ ^ dsi 

a/i^i n y '-'—rx a/ s— sx 
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e). 

_ _ 2 _ 

dri 

/ 8 1 8 \ /''S ds 

itMo , 

-4(1 >J r-ri 

\ 8ri 8s JJ-y.^ J s-si 


2 

dn 

/ 8 ^ ^ V' dsx ^ 


«Mo< 

Jr-^1 

\8ri 8s /'Vi. 3 ^(l_jyiQ) s— Si 


(1+Mo) 

ri(l-Mo) 

r dri _L + r (1+Mo) _ 2 

Jr^ V r-ri \ ^sJJq « 


/ 2sx 

(ri-Si)-Mo(ri+Si) 


— arc tan 


I 2s2 

(ri-Sx)-Mo(ri+sx) 


•J S-S] 


where 


= 


s(l+Mo)-Cy>^ 

1-^ 


Most of these integrals can be readily evaluated to give 


1 f ^ ^ = Jl 



1. f ^ 

^<lo 



— f ^ / ^0-^0 + arc tan f^o*^ ^ 

jtMo \ 1 +Mq V Moto+Xo V to-xo / 


8 


it(l+Mo) 


/ tp-xp ^ _4_ 

V Mq^o+^o 


. 2xo-to(l-^) 

arc sin — ^ + 

(1+Mo) to 


arc sin ^~? o Ho(1+Mq) ] 
to(l-Mo) J 


J- arc Sin 

J to(Mo-l) 


(A6a) 

(A6b) 

(A6c) 

(A6d) 
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16 / tp-xp 

jt^(l+Mo) )J Moto+Xo 


I -EF(t,k») - KE (4r,k») + KF (ilf,k‘) 


2 . X . 4 . 2x-t(l-Mo) . 

— — arc sin — + arc sin — ; — ^ + 

ifMo t JtMo (l+Mo)t 


32K 


y 


1— Xq— M pto 


« (1+Mo) V (l-Mo^)(xo+to) «Mo 


2 - arc cm ^ o. 


to(l+ Mo) 


where 
k» = A-k^ 


(A6e) 




(to+Xo) (l+Mo) 
t = arc sin<y Xo+Moto 

4 


Gi = 




c Vg 

"1+Mo 


dsi 


W{ s-si)(r+si) J 


. (r+si )( 1-Mo) [(s-si)( 1+Mo) ~ c^]+ 2si[s(l+Mo) - c/^r(l-Mo)] 

arc sin — ^ 

[s(l+Mo) - cV2 + Si(l-Mo)][r(l-Mo) - Si(l+Mo)] 


} 


Pitching Wing 


A similar analysis 
results for the loading 


of the pitching wing yields the following 
coefficient: 



(A7a) 
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1/ 8 


9 V 9 o /4 3 «Mo 


, 3/2 


( to^o) + 3y(to-^o)(Moto«o) * 
( 1 +Mo) VMoto■^:^^ 


i/SE') 

1 \1oi 


3(MQtQ+XQ) arc "tan 


Mpto+Xp 

tp-Xo 


I 


•Xp)(Mptp+Xp) + - 


1 Mp(l-Mp) /(tp-xp) 


3 (1+Mp)^ V Mpto+^o 


A 


J [l-(Mo'to‘'’^o)^ t^'^o+Xp) - 1 ] + 


(Mptp+xp) 


arc tan 


A~(Mptp+^ _ „ . / ~^o 

y (tp+xp)-i ^^^yMptp+xp 


1 = -L /( 

<1 VW4 " nMp 


0 


(Mptp+Xo) arc _ 
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y ^~( hd^o'*’^o ^ '*’ 


^ ® - 

r Mo [ 

rtMo 

IVl-Wp^ L 


(I^tp+Xp) tanh ^ yi-(Moto+Xo) 


y(tp Xp)(tp+Xp) 
_ 


y(tp-xp)(Mptp+xp) + 


Mo( 1-M /(tp-Xpf 


(Motp+Xo) 1^ 


arc tan j - arc tan + 

V ^o~Xp , 


to-Xo 


(A7b) 


(A7c) 


(A7d) 
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Mo 




taiih 


— 1 


(to-Xo)(l-Mo) 

(to+Xo)(l+Mo) 


+ ic2 


(A7e) 


where 


-2Mq Vr+Si - 


2s 1 


V r+S] 


— 4si(r+si) + (r(i+Si)[r(l-Mo) + Si(l+Mo) + 2sx(l^4o)] n' 
arc tan - " ■■ + g 

2 72si(r+si)(r-r^)[r^(lHyiQ)-Si(l+Mo)] 


4Mo /n/ r— arc tan 


2s 


1 + 2Mort 


r^(l-Mo)^i(l+Mo) 


^r(l-Mo)— Si(l+Mo) + 


2 

3 1+Mo 


(2si) 


3/2 


S— 5 X 


6 \/2sx 


_1 /ra.(l-Mo)-Si(l+Mo) 'i 

arc tan — 


-1 


(l-Mo)(r-ra) J 


LIFT AND PITCHING-MOMENT COEFFICIENTS 


The lift and pitching— moment coefficients may be obtained by suit- 
able integrations of equations (a 6) and (A7) and are given in the time 
intervals indicated in sketch (aa) by the following expressions: 


Sinking wing 


0 < to < 


1 

1+Mo 



l-to(l-Mo) 


■ Mo 


i - % ^ 


(Mo-2) 


(A8a) 


(A9a) 


1 

1+Mo 


< tf 


< 

~ 1-Mo 
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4 

a. «Mo 


/g-toddi^i). ^ 


1+Mo 


2to(l+Mo)-2 




1+Mo 

4-2to(l+Mo) 


1+Mo 


arc sin 


do(l+Mo)-l 

to(l+Mo) 


^ yto(l+Mo)-l + IS-to(l+Mo)] arc 


1-Moto 

2 


where 


^ol 


1+Mo 


— tr 


t.)^o 


(A8b) 


is given by equation (A6e) 


-ma 


, ^ _ _8_ r r to^(^l8Mgf9Mo^) 2to(Mo-l) 2 ' 

^Mq i5 1+Mq ( 1+Mq ) ^ 


arc tan 


/ 2-to(l-Mo^r r -9to^(l-Mo)^ to(Mo^-4Mo+3) ^ Mo^-f2Mo-3 -j 

y2[to(l+Mo)-l] L 16 2(1 +Mo) 2(1+Mo)2 I 


arc tan 


1 ^ y2[2-to(l-Mo^)][to(l+Mo)-l] 

A/to(l+Mo)-l 1+Mo 


to(l-6Mo+9Mo^) 

16(1+Mo) 


l->-9Mo 

4(1+Mo)^ 


+ ^ yto(l+Mo)-l 

1+Mo 


9to(l-^o) _ 3+3 Mq 
16 '8(1+Mo) 


^ o (xo+Moto)r^^ dxo 
iTmo ~ 


0 < to< 


Pitching wing 
1 

1+Mo 

- ^ 
■ Mo 


^olMo- 


Mo" 


+ to(Mo-l) + 1 


(A9b) 


(AlOa) 
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c 

_ 1 _< 

1+Mo 

c 


where 

ClUj’ 


3Mo 


tn ^ 


{ to^[| (l-Mo)'" + |Mo3 + ^to^(l-Mo)^-| to(l-Mo) + 1 j 


(Alla) 


1-Mo 

8 


It Mo 




‘bo+'to^”l 


tQ(^Mog+2Mo-3) ^Mq^+7Mq+3 
12(l+Mo) 3(1 +Mo)^J 


7toMi^fcto)= (f-Sfa) + 


,(-Mo2+2Mo+ 1) ^ (1-Mo) 


““ "tr 


+ 1 


arc 


tan / ^ — [ - 

V to+toMo"”l \ 

1 Mq 


- - 

4 2 


arc tan 


/ tp+l— tpMp 

tp+toMo-l - 


tanJi 


l-MoV to+toMo-l \ 

l+Mp/Vto+l-toMo 


2 V 

Gs is defined tmder equation (A7e). 


). 


,1-Mptc 


G 2 


-tn+: 


-f- dxp 
2 c 


1+Mo 


(AlOb) 


8 

It Mo 


Mo 


, , tanhVri:^y^i^^^4:i) - 

8 /l=fp y V 1+Mo/ \to( 1-Mo) +1/ 


- "flj; (Mo+1)(-Mo^+4Mo+ 1) - (l-Mp)^ + ^ (l-Mp) - 

^ 

'' *“/ to(l+Mo)-l “ ( 




,, ,.. > , f-to=[(l-Mof +8M0I , to( 9 Mo^Mo-ll) 

yte(l+Mo)-l I ; - 7g(i.Mo) " 


39Mo^+^4Mq+19 


§ 53^1 


to^(2+Mo^) + toMo - 8 

3 S 


Wo 


-to+ 


2° (^+Moto)(^^JdXo 


(Allb) 


1+Mo 
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. APPENDIX B 

EETEEMINATION OF SUPERSONIC, INDICIA!, SECTION LIFT AND 
PITCHING MOMENT CURVES 
THE LOAD DISTRIBUTION 


In the case of the unsteady supersonic wing the expression for the 
velocity potential may he readily obtained hy placing the values of w^ 
given by eq.mtions (5) and (6) in the equation 



vudxi 

*/ (t-ti)^-(x-^i)^ 


(Bl) 


where t is the area on the plan form included in the Mach forecone. 
The loading may then be calculated from the relationship given in 
equation (4). 


Sinking Wing 


The load distribution over the regions A, B, and C shown in 
sketch (aa) axe given by the following expressions: 



1/^^ ^ = Jl 

“V<lo4 Mq 


(B2a) 



Hia 
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4 1 MoXo'+'to . Mo^~l _ _ f Xq 

,— arc cos ^ ^ arc cos I — 

x+Moto nMo \ to 


^ j (B2b) 


1 


Mo2-l 


(B2c) 


Pitching Wing 


1 / Ap A 

For the case of the pitching wing the values of q i ^ j regions 


regions A, B, and C are, respectively. 


5 ( ‘k 


(B3a) 


arc cos ^2^^° + Xp-^Mo 
v^Mq^I ^0+^*0 Mo 


t 2 
^o ^o 


1 / 

<i\ /I L 


ywo^-l 


(xo+Moto) 


(B3h) 


(B3c) 


LIFT AHD PITCHIWG-MOMEm’ CQEFFICIEmB 


The lift and pitching-monent coefficients may be obtained by suitable 
integrations of equations (B2) and (B3) and are given in the time 
intervals indicated in sketch (aa) by the following expressions: 

Sinking wing 


0 < t^ < 


Mo 


(B4a) 


(2-to^) 


(B5a) 



MCA TN 2k03 


—L.<t <-L- 

1+Mo - ° - Mo-1 


^ 7 lir ^ ■ arc cos (to+Mo-'toMo^) + 

, “ L ^ to 




- - r — 

« [ Mo 


T to^ Mo"to“”^ I ^ 

1 =- 1 arc cos . ^ — + 

2 / to 


arc cos 


(to+Mo— toMo^) ■*■ ^ ^ ^/ 1 o^“{ l~Moto ) ^ 


^ to ^ « 




c ’ = — 

% 


Mo^-1 


Pitching wing 


0 < to < 


\ 




c » = / 1 + 

% 3Mo V 2 


Moto^^ 


<t <-i- 

,+Mo ° Mo-1 


. . = 2 r / to! 

'Zq It 1_ Mo i 2 


1 +t^ ) arc cos Msr± + 

2 y to yM“ 


Mo^l 


arc cos (to+Mo-toMo^) + ytoMl-Moto)^ 
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3« [Mo 



Moto^ \ 

2 J 


arc cos 


Mpto-l ^ 1 


arc cos (to+Mo— toMo^) + 


(8— MQto~2to^ — 
6Mo 


yto=-{l-Moto) 



(B7b) 


00 


2 

yMo^-1 


(B6c) 


(B7c) 
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(a) Lift on a sinking wing. 

Figure L- Variation of two-dimensional indicial lift and 
pitching— moment coefficients with chord lengths 
traveled for several Mach numbers. 
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(b) Pitching moment (about leading 
edge) on a sinking wing. 


Figure t. - Continued. 
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Vj’/c / Chord lengths traveled 
(c) Lift on a wing pitching about its leading edge. 

Figure /. - Continued. 




Asymptotic values 






M^^l. 
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VJ/c , Chord lengths traveled 

(d) Pitching moment (about leading edge ) on a 
wing pitching about its leading edge. 


Figure Concluded. 











0 20 40 60 80 too 


Percent chord 
(a) Sinking wing. 

Figure 2. -Variation of two-dimensional load distribution 
with percent chord for a Mach number equal to 
0 . 8 . 
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0 20 40 60 80 too 


Percent chord 

(b)Wing pitching about its ieading edge. 
Figure 2.- Concluded. 
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I I 1 I I ^ I I I I 

0 20 40 60 80 too 


Percent chord 
(a) Sinking wing. 

Figure 3. - Variation of two - dimensional indicial toad 
distribution with percent chord for a Mach num- 
ber equal to i. 2. 
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Figure 3. - Concluded. 
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